Valley-based noise-resistant quantum computation using Si quantum dots 



o 









cd 



O 

(N 
> 

m 

o 
o 
o 

• 

o 



'i 



Dimitrie Culcer/ A. L. Saraiva,^ Belita Koiller,^ Xuedong Hu,^ and S. Das Sarma^ 

^ICQD, Hefei National Laboratory for Physical Sciences at the Microscale, 

University of Science and Technology of China, Hefei, Anhui 230026, China 

^Instituto de Fisica, Universidade Federal do Rio de Janeiro, 

Caixa Postal 68528, 21941-972 Rio de Janeiro, Brazil 

^Department of Physics, University at Buffalo, SUNY, Buffalo, NY 1^260-1500 

^Condensed Matter Theory Center, Department of Physics, 

University of Maryland, College Park MD20742-4111 

We devise a platform for noise-resistant quantum computing using the valley degree of freedom 
of Si quantum dots. The qubit is encoded in two polarized (1,1) spin-triplet states with different 
valley compositions in a double quantum dot, with a Zeeman field enabling unambiguous initial- 
ization. A top gate gives a difference in the valley splitting between the dots, allowing controllable 
interdot tunneling between opposite valley eigenstates, which enables one-qubit rotations. Two- 
qubit operations rely on a stripline resonator, and readout on charge sensing. Sensitivity to charge 
and spin fluctuations is determined by intervalley processes and is greatly reduced as compared to 
conventional spin and charge qubit s. We describe a valley echo for further noise suppression. 



Quantum computing (QC) requires accurate control 
of the two states in a quantum bit (qubit), which en- 
tails long coherence times and ability to rotate a sin- 
gle qubit and entangle neighboring qubits. Solid state 
spin systems are an obvious choice for scalable QC, with 
Si known for its outstanding spin coherence times [iHi], 
thanks to lack of piezoelectric electron-phonon coupling 
[5], weak spin-orbit interaction [6, 7 and most impor- 
tantly, nuclear-spin free isotopes, enabling isotopic pu- 
rification to remove hyperfine coupling [8 . Recently Si 
quantum dots (QDs) have been at the forefront of QC 
research, with theories of spin relaxation in Si QDs fur- 
ther justifying this choice [9l[l0]. Experimental progress 
in Si QC has been made using QDs in Si/Si02 [1THT4] , 
Si/SiGe [151 [IS] ? and donor-based architectures [TTHSS] . 

The two main obstacles to Si QC are dephasing due 
to charge noise and the valley degree of freedom. In lay- 
ered structures only the two low-lying valleys perpen- 
dicular to the interface are relevant, and the interface 
potential gives a valley-orbit coupling A, which deter- 
mines the two-spin exchange [24 . The study of Si valley 
physics has intensified, in experiment [25.-.3Qj and theory 
IIl[3TH35l[37l. 

Here we propose a new QC platform in Si QDs us- 
ing the valley degree of freedom to overcome noise. It 
exploits the closeness in energy of two polarized spin- 
triplet states with different valley compositions. It em- 
ploys available technology in Si QD in MOSFET and 
donor architectures, offering all-electrical control and a 
fundamental perspective on valley manipulation [53j . 

We take an inversion layer of thickness L, with a top 
gate setting the potential at the top Si/barrier interface 
to Vt and a back gate inducing Vb to the back interface. 
The electric field inside the inversion layer is F = (Vb — 
Vt)/L. In a single QD defined in this environment, an 



electron experiences the potential 
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(1) 
with the origin at the top interface. The dot is located 

at {xd^ 0, 0), with a radius a, tti* is the Si in-plane effec- 
tive mass and Uq 6{z) the interface potential with 6{z) the 
Heaviside function. The effective-mass ground state wave 
function is D^{x,y^z) = (j)i:){x^y)ilj{z)u^{r) e^^^^ . The 
lateral envelope function (pnix^y) is a 2D Gaussian, and 
iIj{z) is the solution to the triangular well imposed by the 
interface and the potential difference in the z direction. 
The Bloch functions are u^{r) e*^^^, where the valley in- 
dex ^ = {z^z}^ and k^ = ±/co, with ko = 0.85(27r/a5'i), 
and asi = 5.43A the Si lattice constant. In the basis 
{D^} the single-dot Hamiltonian reads 
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with confinement energy Ed and valley-orbit coupling 

AD = {D,\UoO{z)^eFz\D,), (3) 

The top and back gates independently control the di- 
agonal confinement energy (electron density) and the off- 
diagonal valley-orbit coupling (vertical field.) Gate cross- 
talk can be circumvented using gate compensation tech- 
niques [22 . Experimentally, the lever arm of the top and 
back gates must be obtained in order to determine how 
Vt and Vb control €d and Ajj- The functional form of 
^^(Vt, Vb) and A]j{Vt^ Vb) may be obtained from the- 
ory pi. The quasi-triangular potential in the z direction 
is not exactly solvable. Nevertheless, it is a good ap- 
proximation to adopt the exact expression for an infinite 
barrier E^ = (nV2m^(1.138057reF)2)V3^ [Jg] xhe total 
energy is cd = h^ /m*o? -^ E^ -^ eVr- With E^ fixed, Vt 
shifts the diagonal terms in Eq. ([2|. Such control could 
also be realized in the setup of Ref. 30, where AI-AI2O3 
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Figure 1: Color plot of \Ad\ as a function of the gate volt- 
ages. The dashed lines refer to paths of constant en- If the 
qubit is manipulated along these lines, the back gate voltage 
compensates the change in orbital energy induced by the top 
gate voltage. We particularize our results to Vr > Vb, so that 
the electron is at the top interface at all times. 



multi-gate stacks provide excellent tunability of the ver- 
tical electrical field and inter-dot coupling. 



The valley-orbit coupling A^:) 
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with energies e±. The dominant contribution to /S.d 
comes from Uq. The direct contribution from F^ is 
smaller [1 , yet Fd controls /S.d by manipulating the 
electronic density at the interface, and |Ai)| can be in- 
creased by an order of magnitude. For Si/Si02 barriers 
|AdI = [(27.4±0.02)F] meV, with F in V/nm [38]. Fig- 
ure fl] shows the way Vt^Vb control \Aij\^eD separately. 
We focus on a double quantum dot (DQD), with the 
left dot located at xl = —xq and the right dot at 
xr = xq. The valley eigenstates \D±) are identical in 
both dots if the interface is sharp along the growth direc- 
tion and flat perpendicular to it, or if interface roughness 
is correlated over distances much shorter than the size of 
the QD [2 . We assume the top and back gates can be 
adjusted independently for the L, R dots, so the elec- 
tric field F has different values for the two dots, Fl and 
Fr respectively, while their electrochemical potential is 
constant. The DQD confinement potential is 
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Neglecting interdot couplings due to the interface poten- 
tial (of order IAd/sd^ with / = {L^\R^) <C 1), Ad is 
a constant across the two dots and we diagonalize Hjj 
separately for D = L^R^ obtaining A/,, A^^ and their 



corresponding valley eigenstates. These are orthogonal- 
ized (see Supplement), yielding the states \D±). We use 
only \D±) in the remainder of this work. Initially we take 
Al = Ar^ thus \D±) are identical for L,i?, and interdot 
tunneling preverves the valley eigen-index (i.e. between 
\L±) and \R±), but not between \L±) and |^^).) 

We seek a two-dimensional subspace of the full DQD 
Hilbert space that consists of two states with different 
valley compositions that can be coupled controllably. 
The full Hilbert space includes spin and valley degrees of 
freedom. A high enough magnetic field (~ IT) lifts the 
spin degeneracy, and only the lowest-energy H polarized 
triplet states are energetically accessible, thus below it is 
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The hopping integral i = {L(^\Hq\R^) + (L^i^lVeeli^-R?), 
the exchange integral j 

{LfR^^l\V,,0^\Rf), and 5 = (sl - sr) - 
{u — k) is an effective two-particle detuning be- 
tween dots, with w = {D^'^ Df^\Vee0l^ Df^) = 



and k 
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the on-site Coulomb interaction 
Li ^ R): the direct Coulomb in- 



teraction between electrons on the two dots. In the basis 
fLR^ fRR^ ^-^^ effective Hamiltonian is 
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with \A\tot = \Al\ + \Ar\ and lA^ = \Al\ - \Ar\. In 
the absence of interdot tunneling between opposite valley 
eigenstates, the states T^^ and Tff are decoupled from 
the other three states, thus we do not include them in 
the rest of our analysis. The relevant eigenstates are 



T<_ = {e</^s<' + 2P) [{i/s<) (T^^ - Tff ) + ff^] 

e< = (1/2) [-~s+j - ^Jis+jr+^n 

(8) 
The qubit is defined by the states T^_ and T^^ . Re- 
stricting our attention to the qubit subspace 
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Here hwz = —Sq — j = —Ehmi the Hund-Mulliken en- 
ergy that in one- valley systems denotes the singlet-triplet 
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Figure 2: Two-electron levels in a Si QD. For ^z > 2A, T^^_ 
becomes the ground state and is initialized unambiguously. 



energy difference. In multivalley systems, due to the ad- 
ditional degree of freedom, E^m is not necessarily related 
to the spin part of the wave functions [40 . 

Referring to Fig. [2| two electrons can be unambigu- 
ously initialized into T^^ when J is raised and £r <^ El- 
Then 8 is swept to the (1,1) regime, with T^^ evolving 
into T^_. This setup allows readout by charge sensing 
using a quantum point contact [42]. Raising 5 again, the 
two qubit states have different occupancies, realizing a 
valley blockade. Readout of the charge state of one QD 
(e.g. R) automatically yields the qubit state. 

A workable quantum computer requires tunable x- and 
z- rotations, i.e. the (Jx and a^ gates. A gate electric field 
II z generates the (Jx gate. In the far detuned regime 
the states T^J^^^ and T^^'^ are effectively degenerate. 
Ramping up Fl induces a small difference A^; in the 
valley splitting on the left dot. With an electric field 
difference of 10~^ V/nm, A^; ^ 27.4/ieV, and the gating 
time for a 27r rotation around the X axis is 0.15 ns. For a^ 
rotations Ehm can be altered by changing J, providing a 
az gate of arbitrary strength. Thus, with |A£;| ^ Ehm 
or Ehm ^ \^e\^ a da^ or a a^ gate can be obtained. 

The (Jx and a^ gates also allow a valley-echo experi- 
ment. Both ax and a^ rotations can be gated quickly, 
providing a fast two- axis control, in contrast to e.g. 
singlet- triplet qubit s, where one axis is the fixed Over- 
hauser field gradient. With the available fast gates, ad- 
vanced dynamical decoupling techniques [43l |44] can re- 
store coherence. If A^; < Ehmi ^ rotations may be 
implemented through valley resonance (cf. Ref. Tl) An 
alternating Asit) at a frequency resonant with Uz en- 
ables Rabi-like flips of the valley qubit. 

Initialization requires only a magnetic field, while in 
the far detuned regime, where operations are performed, 
the magnitude of A is irrelevant. Coherent rotations can 
be demonstrated without precise knowledge of A, even 



for A < ksT. Knowledge of A is required for controlled 
rotations, whose characteristic time scale and coherence 
time may be estimated by Ramsey interferometry. 

The universal set of quantum gates is completed us- 
ing a superconducting transmission line resonator to pro- 
vide entanglement. Superconducting resonators have ex- 
ceptionally small mode volumes for strong couplings to 
qubits, and very high quality factors. They can mediate 
interactions between singlet-triplet qubits [45l [46] and 
can couple valley qubits. An array of DQDs can each 
be coupled to a superconducting transmission line res- 
onator, with only the right dot coupled out of each DQD. 
The (quantized) resonator voltage is an addition to the 
detuning, and the effective interaction is 



Hres=g{a + a^)\T^!!){T^ 



^nRR\ 



(10) 



where a is the annihilation operator for the lowest energy 
mode of the resonator. Hres is mapped onto a Jaynes- 
Cummings Hamiltonian. For Si shallow donor electrons 
coupling to Si cavities, g has been estimated at 30 MHz 
J47] , and for QDs we expect an improvement of an order 
of magnitude. Type-II superconducting alloy striplines 
have high Bc2'- NbTi, which is already used, has 5c2=15 
T, and can be operated at the magnetic fields required 
here to freeze out the spin degree of freedom. 

Valley qubits have good coherence properties. Firstly, 
elaxation mediated by a generic phonon potential Vphn 
is quantified by Ti, and strongly suppressed, because 



1 



27r 



^ \{D^\V,Hn\D-)\'S{e^ - s_ - M (11) 



requires phonons with energies hu ^ 2A and k of the 
order of ko. Optical phonons bridge the difference in k 
but have energies ^ 2 A, whereas acoustic phonons with 
energies ~ 2 A have wave vectors <C ko [48]. For the same 
reason, the qubit does not relax to the lower- lying T^^ 
states with orthogonal valley composition. 

Valley qubits are immune to charge noise due to dan- 
gling bonds (DBs), which hamper conventional charge 
and spin qubits [49]. Since a large wave vector 2ko sep- 
arates the z and z valleys, intervalley matrix elements 
are only appreciable for real-space interactions sharp in 
the z-direction. For interactions non-singular in z, ma- 
trix elements between \Dz) and \Dz)^ likewise between 
\D-\-) and |^-), are suppressed. The Coulomb potential 
v{x^y) of a DB in the QD plane has intervalley matrix 
element (v) = {Dz\v{x^y)\Dz) = 0, since the overlap of 
valley states vanishes (see Supplement.) For a DB not in 
the plane, the intervalley matrix element of v(x, y) is still 
suppressed, since quasi-2D screened Coulomb potentials 
remain long range. Thus the ax gate is immune to noise. 

Noise in the a^ gate could lead to pure dephas- 
ing [49 . Variations in Ehm are given by AEhm = 
{dEHM/di)Ai ^ {dEHM/dS)AS. In the far-detuned 
regime, where the qubit is operated, there is little vari- 
ation in t [49]. With the ax gate off, Hqubit = {Ehm + 



^Ehm) (^z and consider the decay of the off-diagonal ele- 
ment pi2 of the density matrix. Random telegraph noise 
causes pi2 ex e~^^^ cosr^t, with rj ^ AEhm/^ and r the 



switching time of the fluctuator, which can reach ms 
For 1// noise, pi2 ex e~^^*\ where 
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where S^{uo) is the spectral density of fluctuations in 5 
and uoq a cut-off (the inverse measurement time.) Since 
the two qubit states have the same envelope functions, 
Ehm is effectively independent of 5 (Fig. 2 of Supple- 
ment). The sensitivity of the a^ gate to charge noise can 
be reduced to any desired level by reducing dEHu/dd. 

Charge noise is dominated by the low-frequency part 
of the spectrum, is long wavelength and cannot differ- 
entiate valleys. In conventional charge qubit s with dif- 
ferent charge distributions, charge noise and electron- 
phonon coupling differentiate the states and lead to de- 
phasing. We have devised a charge qubit that is im- 
mune to charge noise, and whose coherence can be further 
improved through valley-echo. The valley-based qubit 
states have the same spin configuration, so this architec- 
ture is also insensitive to magnetic noise. (A forthcoming 
publication [50 obtains an intervalley dephasing time of 
/is, rather than ns for GaAs charge qubits.) These ar- 
guments still hold if the effective mass approximation is 
relaxed by including /c-states in the vicinity of ±/co- 

Charge defects in the dielectric and interface roughness 
affect the vertical field, and may cause |Ai,| ^ \^r\' The 
operation of the qubit is unaltered, since a change in |A| 



still induces mixing between the eigenstates of Ht- Any 
difference |Al| — |Ai^| can be offset by gate control of A 
prior to operating the qubit. The proposal is robust with 
respect to unavoidable differences between the two dots. 

The top gate may alter the phase of A by a small 
amount [1 , inducing tunneling between the qubit space 
and the states T^^. Yet T^^ are separated by 2|A|, and 
the intervalley tunneling matrix element is a fraction of 
the intravalley one (which is ~ tens of jieV [40 , an order 
of magnitude less than |A| as measured recently [ 28H3Q] .) 

We have devised a new valley-based QC platform in Si 
QDs that offers full electrical control with current tech- 
nology. The qubit is initialized using a Zeeman field, 
rotated using gate electric fields, and read out via charge 
sensing. Coupling to a superconducting resonator allows 
entanglement. The platform is tailored to Si, with the 
valleys perpendicular to the interface, and A controlled 
by a top gate [l] [27] . Valley-based qubits in carbon QDs 
must be studied separately [5TJ|52]. Charge-based QC 
schemes, being easier to control, yield insight into qubit 
manipulation and coherence, and spur better design. 
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SUPPLEMENTARY MATERIAL: ONE AND 
TWO-ELECTRON STATES OF A 
MULTI- VALLEY QUANTUM DOT 



Figure 3: Spin-polarized ^^ triplet levels in a DQD, as deter- 
mined in Ref. 2, with a valley splitting of O.lmeV. 



Each dot is located at Rd^ with D = L^R. The 
ground-state wave functions D^ satisfy {T -\- Vd) D^ = 
£d D^^ with T the kinetic energy operator including the 
effective mass anisotropy (m* ^ rriz^ with tti* the in- 
plane effective mass as in the text), 



ikpz 



D^ = <PD{x,y)^{z)u^{r)e''^ 



(13) 



The lattice-periodic function u^{r) = "^Zk ^k^^^'^ •, with 
K the reciprocal lattice vectors. The envelopes (j)D{x^y) 
are Fock-Darwin states 



(t^Dix.y) 
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The EM A equation satisfied by '^(^) has been given in 
Refs. [11 O It can be solved numerically or analytically. 
For analytical calculations, a convenient choice is the 
variational wave function [3], 



iIj{z) = Me^.z < 



(15) 



N {z^ zo)e ^\z > 0, 



where k^ = J ^^2 ° , and continuity of the wave function 
at the interface requires M = Nzq. The only variational 
parameter is ksi^ with the others given in Ref. |2| The 
overlap of states from different valleys is 
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Since Iz^z <^ 1 by many orders of magnitude we take 
Dz, Dz to be orthogonal. The overlap is zero in bulk Si, 
and it remains zero if the envelopes are wide and smooth 
compared to the lattice constant. This is the case in our 



numerical analysis. The appearance of Iz^z is an artifact 
of the trial analytical wave function we have suggested. 
The overlap / = {L^\R^) is not zero. It is therefore con- 
venient to construct single-dot wave- functions that are 
orthogonal, as was done in Ref. |2l as Z^ = . ^~^ ^ = 



and R<: 



Rj-gL^ 



, where g = (1 — y/l — P)/l^ so that 



H — I ^ 

(,Ri\li) =0. We define 

A= {b^\V^Q{z)^eFz\b_^) (17) 

t = {L^\Ho\R^) + {L^L^\Vee\L^R^) = h + s. 

We subsequently orthogonalize L± and R± obtaining 
Z± = L^-9R± and R± = f^~^^^ . These are the 

^1-2/^+5^2 ^ yi-2^^+^2 

states used throughout this work. 

Two electrons are initialized on one dot. The four 
lowest-energy two-particle spatial wave functions are, 

qRR _ 5(1) 5(2) 

S^^ = {l/V2){R^l^R^^^^RfR^l^) (18) 

Tf^ = {l/V2){R^l^R^^^ - RfR^l^), 

with the superscript i denoting particle i. The following 
two-electron matrix elements are needed 

i = 0i^RfK.\LfRf) = 0i^R^^l\V^e\L^^lRf) 

u = {bfbf\v.Abfbf) = {bf^b^^lw^Abf^b^^l). 

(19) 




Figure 4: Bloch sphere representation of the qubit states T^_ 
and T^^ . The Gx gate is provided by the top-gate induced 
difference A^; in the valley-orbit coupling between the two 
dots. The Gz gate (not shown) is provided by the detuning. 



In the basis {S^^, 5f ^, Tf ^, ^f f } the Hamiltonian is 



H 



Idot 



2sn 



/- 
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ml 



(20) 



A Zeeman field can separate the polarized triplets, and 
for a Zeeman energy Ez > 2 A, the lowest level is the 

polarized triplet T^ . 

When both dots are considered, we will only re- 
quire two-electron spin- triplet states. For an asymmetric 
DQD, in which the potential of the left dot is raised, we 
may neglect the LL states. The remaining triplets are 



pLR 

^-±± 

pLR 



(l/v/2) {L^l^R^i^ 
(1) 5(2) 



{i/V2){L^l^R^^^-L^i^R 



^(2)5(1)^ 



(21) 



ff^ = (1/V^) (^V'^^-^ - ^v^^' 



?(2)5(i)^ 



In this work we only consider the polarized H triplets, 
thus in all of the above it is understood that e.g. T^^ = 
^±± ' ^^^ ^o forth. The effective Hamiltonian in the 



basis {f^^,ff^,ff^,f^f,ff^} is 



Hj 



/-\A\tot-j 0\ 

|A|,,,-J 

\A\e -j i 

-j -\A\e -i 

i -i -5j 



(22) 



When the Hamiltonian Ht of Eq. ( 22 ) is diagonalized for 
L 



0, the eigenstates are T^^ an3~ 



T>_ = 



f<_ 



fiVlf^, 



,>2 



2*2 



rpanti 
f2 V ^0 



-0 



(W2f^^ 



rpanti 
-<2 + 2*2 V 4 



_fRR 



niifi 



(23) 



rpSym -*- /rpLR , rpLR\ 

\f2 



rpanti 



1 



V2 
The energies are 



(n 



LR 



Tt'^). 



-6 + j 


+ ^i} + ~jY 


+ 8*2 




2 




-5 + j 


-^{5 + ~JY 


+ 8*2 



(24) 



Figure [3] shows the five energy levels as a function of the 
dimensionless detuning in the case of no valley coupling 
difference across the two dots. The qubit is defined by the 
ground and second excited states on the right of Fig. [3) 
with energies e^ and 0. Figure [i] illustrates the qubit 
explicitly on the Bloch sphere. 

Valley-echo removes the effects of inhomogeneous 
broadening and the associated dephasing. The state T^^ 

is initialized and swept to T^_ , then A^; is set to a fixed 
value, controlling the cFx rotation. After a specified time 
r the detuning is swept so as to apply a fast 7r-pulse of 
the (Jz gate. Readout is performed at time 2r, when the 
qubit returns to its initial state. 
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